Non-Markovian theory for the waiting time distributions of single electron transfers 
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We derive a non-Markovian theory for waiting time distributions of consecutive single electron 
transfer events. The presented microscopic Pauli rate equation formalism couples the open electrodes 
to the many-body system, allowing to take finite bias and temperature into consideration. Numerical 
results reveal transient oscillations of distinct system frequencies due to memory in the waiting time 
distributions. Memory effects can be approximated by an expansion in non-Markovian corrections. 
This method is employed to calculate memory landscapes displaying preservation of memory over 
multiple consecutive electron transfers. 



I. INTRODUCTION 

Detection of single electron transfers through quan- 
tum systems such as quantum dots has become exper- 
imentally feasible ii 2 -! 3 -^ Theoretical investigations on the 
underlying statistics were mostly obtained in terms of 
higher cumulants, e.g. noise and skewness, by using gen- 
erating function techniquesi^&^&^ i 10 ' 1 ^ 12 ' 13 ' 14 Expan- 
sion of the higher cumulants in non-Markovian correc- 
tions has revealed significant memory effects in quan- 
tum dotsi^ when strong Coulomb interaction, 1 ® phonon 
bath— or initial correlations^ are present. 

Statistics based on waiting time distribution (WTD) 
provides additional information on the system. Higher 
cumulants can be derived from the WTD 19 , but not 
vice versa. Waiting times were recently utilized to 
analyze single electron transfers in the Markovian 
regime, for example, in double quantum dots,— single 
molecules ) 21 ' 22 single particle transport^ and Aharonov- 
Bohm interferometers. 23 Non-Markovian treatment of 
WTD has shown significant features in photon counting 
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statistics. 

Non-Markovian effects are induced by a small bias volt- 
age or a finite bandwidth of the system-electrode cou- 
pling. While the former can be eliminated easily, the 
latter scenario is given by the setup of experiment. In 
order to explore both regimes, a non-Markovian Pauli 
rate equation based on a microscopic description of the 
electrode-system coupling using a Lorentzian spectral 
density is derived. It can be utilized for a variety of 
system, such as single molecule and quantum dots. 

A formal connection of the WTD with the shot noise 
spectrum of electron transports through quantum junc- 
tions has been established in the Markovian regime^ 
The non-Markovian shot noise spectru m 15 i 25 provides a 
more accurate description of the signal and its relation 
to the physics of the junction than a Markovian version, 
since it reveals several distinct intrinsic system frequen- 
cies. 

In this paper we derive a non-Markovian theory for 
the WTD of single particle transfer trajectories based 
on the derivation of a non-Markovian microscopic Pauli 
rate equation. It provides a general framework to 
study non-Markovian electron transport through many- 
body systems and allows us to distinguish between non- 




u 
■a 

o 

H 
O 

a 
"qj 

(50 



FIG. 1: Illustrated set-up of the DQD in series and notation 
of the important parameters. The charge state of the DQD 
is measured by the quantum point contact (QPC) which pro- 
vides an electron transfer trajectory. 



Markovian effects due to intrinsic properties of the sys- 
tem, finite electrode-system coupling band-width and 
small bias voltage. The WTD is evaluated in time do- 
main by perturbation theory leading to non-Markovian 
corrections.— We shall analyze the effect of memory on 
consecutive electron transfers through double quantum 
junctions (DQD), see Fig.[TJ and demonstrate the influ- 
ence of many-body Coulomb coupling on memory land- 
scapes displaying the memory that is preserved in the sys- 
tem for several consecutive electron transfers. The non- 
Markovian spectrum is obtained from a Laplace trans- 
formation. The results reveal that the non-Markovian 
spectrum of the WTD provides similar information con- 
tent which qualifies it as an alternative method to the 
non-Markovian shot noise spectrum. 

The paper is organized as follows. In section [TlJ we 
present the derivation of the non-Markovian rate equa- 
tion. The expressions for the non-Markovian WTD are 
shown in section IIIII The formalism is applied to the 
DQD system and the results are given in section Hvl We 
conclude with a summary and outlook. 



II. NON-MARKOVIAN RATE THEORY OF 
QUANTUM TRANSPORT 

A. Hamiltonian 

Consider a junction consisting of a DQD in series as 
the system, two electron reservoirs, and the respective 
system-reservoir coupling, as shown in Fig.[TJ The total 
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Hamiltonian assumes Ht — Hs + Hr + H$r- The system 
part describes the DQD which is modeled by 



and 
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E s n s + Uhifi2 — A(c[c2 + c\c\ 



(1) 



Here, n s — c\c s is the electron number operator of quan- 
tum dot s = 1 or 2 with orbital energy E s , and U specifies 
the Coulomb interaction between two dots. The reser- 
voirs of left and right (a = I and r) electrodes are de- 
scribed by 
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The system-reservoirs coupling responsible for electron 
transfer between the system and the electrodes is 
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The electron creation (annihilation) operators c\ (c s ) 
and c^ aq (c aq ) involved in Eqs. (HJ)-© satisfy the anti- 
commutator relations. In this system, single electron 
transfer trajectories can be obtained from the charge 
state of the DQD that is constantly measured by a quan- 
tum point contact (QPC). Such a configuration was em- 
ployed in experiment 4 operated at a small bias voltage. 



B. Generalized non-Markovian rate equation 

We now turn to the non-Markovian rate equation. Let 
p(t) = tr RpT(t) be the reduced system density opera- 
tor. The total density operator is assumed to be ini- 
tially factorisable into a system and a reservoir part, 
Pr(to) = p(to)pR(to), and the system-electrode cou- 
plings are assumed to be weak. Using the standard ap- 
proach, one can readily derive a non-Markovian quantum 
master equationi 26 i 27 ' 28 For the present study, we adopt 



T ( S q ] = TF'T$*> for simplification. A rotating wave ap- 
proximation to the cross coupling terms between the sys- 
tem orbitals is not required here since each electrode is 
coupled to one orbital site only^S We denote the sys- 
tem Liouville operator Cs- = [Hg, ■] and set H = 1. The 
resulting quantum master equation in the time-nonlocal 
form reads 20,28 
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The reservoir correlation functions, 
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Ci-\t) = y2\T^\ 2 (c aq (0)cl q (t)) Roi , 
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contain the properties of the electrodes. Here, c4 5 (t) = 
e^'c^e- 1 ^' and (0) Ra = trR a {OpR a }, with p Ra 
being the density operator of the bare electrode a under 
a constant chemical potential p a . Physically, C^t\t) 
describes the process of electron transfer from the a- 
electrode to the system, while d (t) describes the 
reverse process. These two correlation functions are 
not independent; they are related via the fluctuation- 
dissipation theorem. 

Electron counting experiments are operated either in 
the large-bias limit in order to achieve a directional tra- 
jectory of single transfer events 1 "^^ or at small bias in or- 
der to realize transfer against the direction of the bias. 4 
Non-Markovian effects are either due to small bias or 
finite band-width. In order to study both regimes, we 
derive a rate equation by projecting the master equa- 
tion ((4]) into the Fock space of system and by considering 
only the population part p rn = p mm . Some simple alge- 
bra leads from Eq. (T5|) to the non-Markovian Pauli rate 
equation 
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drK mn (t - r)p„(r). 
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Here, u> mn = E m —E n is the transition frequency between 
two Fock states; 



r^HTj">| 2 |<m| C t|n)| 2 , 



(8) 



with s = 1 or 2 for a = I or r, respectively, is the 
state-dependent non-Markovian system-reservoir cou- 
pling strength. As inferred from Eq. (jSJ), Fmn ^= only 
if \m) has one more electron than \n). We can therefore 
identify the rate kernel elements involved in Eq. ([7]) with 
three physically distinct contributions 

K(t) = Y,[K {a+) (t) + K^-Ht)} + K (t). (9) 

a 

K^ a+ \t) and K^ a ~\t) realize an electron transfer in and 
out of the system through the a-electrode, respectively. 
They summarize the off-diagonal matrix elements of the 
transfer rate kernel K{t) in Eq. ([7]), 

K ( «+\t) = rWc< + )(< - r ) e -<W + cx . ; (10) 
K^-\t) = r<$C<->(t - r)e-*-» 4 + c.c. (11) 
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Ko(i) summarizes the diagonal matrix elements of K(t) 
and leaves the number of electrons in system unchanged. 
These diagonal elements satisfy 

(K ) m = -J2 K'W + Ktn\t)] . (12) 



For the Lorentzian spectral density model [Eq. (|A.2[) ]. 
where the reservoir spectral density assumes the form 
J a {ijj) = 7„/[(w - VL a ) 2 + 7^], we obtain for the off- 
diagonal rate kernel elements the following expressions, 

Kffi(t) = 2T%1 {e-^[ai cos(fi° „*) - 6± sin(^ m t)] 



fe=l 



(13) 



Here, Wk = (2fc — 1)tt/ (3 is the fermionic Matsubara fre- 
quency, while f^„ = fl a - u mn and /i° „ = [i a - ui mn . 
The coefficients a^, b^, c^, and are all real, given 
explicitly in Appcndix n by Eq. (|A.10[) . The first term in 
the curly brackets of Eq. (fl"3j) reflects the spectral prop- 
erties of the electrode-system coupling, while the second 
term arises from the decomposition into Matsubara fre- 
quencies which induces memory effects due to small bias 
voltages. From the expressions one can infer that large 
7 Q , wide bands, and large Wk, high bias, cause a fast de- 
cay of the transfer rates in time. The decay is responsible 
for the memory loss in the system. 

The non-Markovian Pauli rate equation ([7]), in terms 
of the population vector p(t) = {p m (t)} and the involved 
transfer matrices, is 



p(t) = [ drK(t - t)p(t). 

J t 



It reads in Laplace frequency domain 

sp(s) -p = K(s)p(s). 
The corresponding electron transfer rates are 



(14) 



(15) 



±y mn \ a ) ^ x rnn\ 



(a) ( a a( s +J a ) -feq»£» 
' (* + 7a) 2 + (fi£») 2 



+E 



k=l 



( S + tu fe ) 2 + (,i«J 2 



-}. (16) 



The derived non-Markovian rate equation formalism is 
based on a microscopic description of the electrode- 
system coupling, and is valid for arbitrary bias and tem- 
perature. Compared to the quantum master equation 
in the same regim o 27 ' 28 , the exclusion of the coherence 
makes it numerically feasible to calculate multilevel sys- 
tems such as large molecules j^ 2 - This allows to include 
non-Markovian effects in large many-body systems, e.g. 
quantum-chemistry calculations, since the properties of 
the molecular-junction enter only through the couplings 



rim and the fitting parameters of the Lorentzian spec- 
trum. 

To rate equation (|14[) . the Born-Markov approximation 
can be applied by separating the integration variables 
and extending the upper limit to infinity in Eq. (|14j) . The 
resulting integration over time, 

/>oo 

W^ = dtK^(t) = K^(s)\ s=0 , (17) 



gives the Markovian electron transfer rates. The second 
identity is via the Laplace domain rate equation (|15p , by 
which the Born-Markov approximation amounts to the 
zero frequency contribution. 



III. NON-MARKOVIAN WAITING TIME 
DISTRIBUTION 

A. Statistics analysis 

We consider two consecutive electron transfers con- 
tained in a time series as illustrated in Fig.[TJ An electron 
entered the system from the left electrode at an earlier 
time to is detected at time t leaving the system through 
the right electrode. No other electron transfers are de- 
tected in between. The joint-probability for the consec- 
utive electron transfer events is 



P(t) = ((W^G{tM)W (l+) p{t ))). 



(18) 



Here, ((■ • •)) denotes the sum over the final system states. 
We assume that the transfer events are instantaneous 
compared to the time-scale of the system propagation 
in between as shown in Fig.[T] Therefore we have used 
the Markovian forms of rate matrices, for the ascribed 
two consecutive events. This assumption is reasonable 
in accordance with electron counting experiments, where 
typical waiting times are long compared to the fast trans- 
fer events The memory of the system is contained 
in G(i, to), the non-Markovian propagator of the system 
from to to t in absence of transfers. It is therefore as- 
sociated with the diagonal rate matrix Kq of Eq. (|12[) . 
satisfying 



dt 



G(Mo) 



dr'K (t-r')G{T',t ). (19) 



For the given two-event case, the joint-probability is 
equivalent to a waiting time distribution. 2 ^ 22 

Now consider the event of an electron transferred 
into the system and the subsequent waiting time be- 
fore any other transfer takes place, also referred to as 
survival probability. In the present notation it is given 
by ((G(t,t )W^p(t ))). While the joint probability is 
subject to the nature of the second transfer, the spe- 
cific form of the second event is irrelevant to the survival 
probability. Consequently, we introduce the survival time 
operator 



Z^ a± \t,t ) =G(t,t )W( a 



(20) 
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If memory is absent, the survival probability is indifferent 
from the previous waiting times. To study the memory of 
a previous survival time that carries on into the follow- 
ing survival time, we introduce two-time joint survival 
probabilities of the form 

QCT-a.T-i) = ((Z('-)(T 2 ,r 1 )^+)(r 1 ,t )p(i )}}. (21) 



B. Non-Markovian corrections 

The formal solution to the propagator in Laplace do- 
main is given by 



G(s) 



1 



s-Ko(s) 



(22) 



The complex Laplace frequency s — 7 + ioj is associ- 
ated with the system residing in its state. The bilateral 
Laplace transformation reduces to a Fourier transforma- 
tion by setting 7 = 0. Since K (s) is strictly diagonal in 
the many-body eigenspace of the system, the matrix in- 
version required in Eq. (|22[) can be efficiently carried out 
for large systems. 

The technique of expanding the propagation into non- 
Markovian corrections has been applied to electron trans- 
port recently ! 16 ' 17 ' 18 Here we apply it to the WTD. Let 
us first express Eq. (|2"2"|) by its series 



n=0 



[Ko(s)Y 



(23) 



Assuming the derivative d 7 s n [Ko(s)] exists for all m, 
the kernel can then be expanded into a Taylor series 
[Ko(s)] n = E™=o Wo(*)H s =o^. Thus, 



G(s) 
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d™[K (s)F\ s = 

ml v n+l—m 



(24) 



Now we apply the inverse Laplace transform x{t) = 
i n ™ ^ st \ ^ Q sw jtch back into time domain. 

27TI Jj — lOO V ' 

One can simplify the poles by using m — n which ne- 
glects the transient terms . 16 ' 17 ' 18 We obtain 



n=0 



]n e K (s)t 



s=0 



n=Q 



G {n) (t), 
(25) 

with GW(i) denoting the individual term involved, 

where G^°\t) = e^ ^* | s =o describes the Markovian dy- 
namics. The first identity of expression (|23|) is asymp- 
totically exact since the dynamics is reduced to the poles 
m = n. The WTD can also be expressed in terms of 
P(t) = Y, n P {n) {t), with P(°'(t) denoting the Markovian 
contribution; so can the survival probabilities. 
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FIG. 2: The relative non-Markovian spectrum F(uo) of the 
WTD. The parameters used for the four panels are given as 
follows. Upper left panel (a): U = 0, AE = 0, V = 1.0, 
A = 1.0,5.0,8.0. Upper right panel (b): U = 0, AE = 0, , 
A = 5.0. V = 1.0,2.0,5.0. Bottom left panel (c): AE = 0, 

V = 1.0, A = 5.0 , U = 1.0,2.0,4.0. Bottom right panel (d): 

V = 1.0, A = 5.0 , U = 0.0, AE = 0.1, 1.0, 2.0. 



IV. DEMONSTRATION AND DISCUSSION 



We employ a non-Markovian rate equation to calcu- 
late the two-electron system as illustrated in Fig.[TJ This 
system resembles the counting experiment conducted in 
Ref.y. Here, the DQD provides a total number of four 
eigenstates: the unoccupied (|0)) two single-occupied 
and |2}), and one double-occupied (|3)), with the en- 
ergies of e = 0, e 1/2 = i(S 1 + J B 2 )T v /i(^i - E 2 f + A 2 , 
and 63 = Ei + E2 + U, respectively. The equilib- 
rium of the chemical potential of the electrodes is set to 
/i cq = (E1+E2) /2. For numerical demonstrations, we use 
the numbers in accordance with recent electron counting 
experiments of electron transfers through quantum dot 
systems at small temperatures? 3 A coupling strength of 
r = 10 4 Hz serves as the unit for all values. This is equiv- 
alent to an energy unit of [E] = 10 4 h = 6.63 x 10~ 30 J, 
and a time unit of [t] = 0.1ms, which is the typi- 
cal time scale of waiting times in quantum dot count- 
ing experiments^ We also use a low temperature of 
T = 2 x 10 4 [i?] = 10 mK. If mentioned, we set a small 
energy detuning of AE — E\—Ei in order to deduce spe- 
cific frequencies of the systems. The bandwidth 7 is set 
sufficiently large in order to neglect the finite bandwidth 
effects; thus the non-Markovian effect is studied in the 
wide band region. In addition, the Lorentzian spectral 
densities are aligned to the orbitals of the system. 
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A. Transients and Fourier spectrum of WTD 

Figure [2] shows the relative non-Markovian spectrum 
of the WTD represented by 



two-times joint probabilities 
)(«) 



F(u) 



\P(co) - P«>Xu)\ 



(26) 



It is noteworthy that F(u>) is independent of the system 
reservoir coupling strength parameter T. The WDT spec- 
trum reveals several frequencies that are present in the 
transient oscillations. These depend only on the internal 
transfer rate A, Coulomb coupling U, and bias voltage 
V. The specific values of the parameters are given in the 
caption of the figure. Figure^a) shows the main charac- 
teristics of F(u>), consisting of two overlapping sub-peaks 
centered around the value of A. Changing the value of 
A leads to the shift of both sub-peaks equally by A. 
In Fig.[5{b), A is kept constant and the bias voltage is 
varied. We find that the splitting of the two sub-peaks 
is determined by the applied voltage. Labeling the two 
peaks with ±, respectively, we can deduce the following 
relation for the corresponding characteristic frequencies. 
lo± = | A ± V/2\. In the presence of Coulomb interac- 
tion, we observe an additional double-peaks feature at 
\U — A ± V/2\, as demonstrated in Fig.^c). This is sim- 
ilar to a non-Markovian shot noise spectrum^ where 
a finite Coulomb interaction U induces also additional 
peaks due to the energy gap between the two-particle oc- 
cupation state and lower states. On the other hand, the 
orbital detuning does not induce additional peaks in the 
double quantum dot in series as shown in Fig.[2fd). 

Oscillations of Rabi frequency, which were observed in 
parallel DQD systems ; 19 ' 20 ! 23 are however absent in the 
present series DQD system. In the parallel cases, the 
transport proceeds via two channels, and the Rabi oscil- 
lations in the WDT can be observed as the consequence of 
quantum mechanical interferences ^ 9 ' 20 ' 23 It is also noted 
that the information contained in the spectrum of the 
non-Markovian WDT is mostly equivalent to a measure- 
ment of the non-Markovian shot noise spectrum. For this 
purpose, the WTD can be considered as an alternative 
approach to the shot noise spectrum measurement. 



B. Memory landscape of consecutive waiting times 

The expansion in non-Markovian corrections, Eq. (|25p . 
can be readily employed to calculate the two propagators 
involved in the two-times joint probabilities defined by 
Eq. 1(211). Denote 



xGW(ri-to)^ ,+ 'p(to)». (27) 

A memory landscape of the system can be calculated 
by the difference between non-Markovian and Markovian 



^'{T 2 ,n) = — w . (28) 

The order n of the perturbative expansion in non- 
Markovian corrections has to be chosen in accordance 
to the parameters in order to assure satisfactory conver- 
gence. We find that the summation to the fourth non- 
Markovian contribution already converges sufficiently for 
the given parameters. As the memory in Eq. (|28|) de- 
cays, the relative non-Markovian landscape lS n ' (t2,ti) 
converges to zero. FigureO shows the memory landscape 
of two survival times related to two consecutive electron 
transfers through the left electrode. It visualizes how 
memory of the waiting time T\ after the first transfer is 
carried over into the waiting time t 2 following the second 
transfer. 

We find that the non-Markovian effects are small for 
the given parameters in case the DQD is coupled sym- 
metrically to the electrodes. This is due to the relatively 
weak coupling of the DQD to the electrodes which is re- 
quired in present counting experiment in order to resolve 
single electron transfers on the measurable timescales. 
It is observed that a stronger coupling to only one elec- 
trode induces significantly larger non-Markovian effects 
as shown in the left panels of Fig. [3] 

In general, fast electron transfers arc necessary in order 
to observe significant non-Markovian effects. The devi- 
ations for Ti, t 2 approaching zero from the Markovian 
value are due to the truncation of the transients in the 
derivation of the expansion. The expansion follows the 
general trend of a numerically exact solution. Both solu- 
tions overlap after the transients have decayed. However 
this causes relatively large inaccuracies for n and t 2 close 
to zero. 

There is an interesting dependency of the non- 
Markovian effects in the memory landscape on the 
Coulomb repulsion U . By comparing the upper panels of 
Fig. [3] where Coulomb repulsion is absent, with the bot- 
tom one, where a large U induces a Coulomb blockade 
regime, we observe that memory decays faster with t 2 in 
the Coulomb blockade regime. This can be explained as 
follows. In the second regime, only a single electron can 
occupy the DQD and the double occupancy state does 
not provide memory for the second survival time lead- 
ing to an overall smaller non-Markovian contribution. In 
this case, only one possible trajectory in the left to right 
direction is possible. An electron enters the unoccupied 
DQD at time T\ = and leaves it at time t 2 = 0. 

The memory is preserved during t± by the single elec- 
tron inside the DQD. However, after the electron has left 
a junction, the memory of its trajectory is lost rapidly 
since no other electron can serve as a messenger inside 
the DQD thus leading to comparatively short survival 
times t 2 where memory is present. In the regime where 
Coulomb repulsion is neglectible, a second electron can 
occupy the junction along the described trajectory, which 
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and observe that the non-Markovian effects increase with 
the coupling strength. This observation is expected to re- 
main true based on general Pauli rate equation dynamics. 
In other words, the non-Markovian effects are mainly vis- 
ible for stronger couplings. 

The employed microscopic non-Markovian rate equa- 
tion provides a general framework to study similar sys- 
tems and allows us to distinguish between non-Markovian 
effects due to intrinsic properties of the system, finite 
electrode-system coupling band-width and small bias 
voltages. It can be combined with quantum chemistry 
calculations that can calculate the employed parameters 
for molecules and their binding to the electronic bands 
of the metal electrodes. The approaches derived for the 
non-Markovian WTD are general and can be applied to 
a variety of processes in physics, chemistry and biology 
that are described by rate equations. 



20 40 60 80 100 20 40 60 80 100 

T 2 [ms] r 2 [ms] Acknowledgments 



FIG. 3: L\ I memory landscape of consecutive survival times. 
The bandwidth is large and a finite bias of V = O.lfctT is 
applied. The left coupling strength is F (I) = 10 4 Hz. The 
upper panels are calculated in absence of Coulomb coupling, 
U = 0, bottom panels display the Coulomb blockade regime, 
U = oo. Left panels are calculated for a symmetric system, 
p(() _ p( r ) j n the right panels, a stronger coupling strength 
is applied to the right electrodes F' 1 "' = W^Hz. 

is represented in the model by the presence of an occu- 
pied double occupancy state. The presence of the second 
electron preserves the memory during ti after the other 
electron has left the junction. 



V. CONCLUSION 

We find that non-Markovian effects are small in the 
regimes of recent single electron counting experiments. 
The sampling rate of current experiments is slow, a re- 
quirement which is imposed by the detection process of 
single electron transfers with currently available technol- 
ogy. This verifies the reason that the Markovian approx- 
imation of previous studies considering FCS or WTD is 
reasonable for the previously investigated systems. 

Non-Markovian effects in the electron transfer statis- 
tics have to be taken into consideration for stronger 
electrode-DQD couplings, which then also requires faster 
sampling rates or a strongly asymmetric system. For ex- 
ample they affect the decay rates of the WTD which are 
directly related to the electronic structure of the system 
in junction^ Non-Markovian effects also induce several 
oscillations with characteristic system frequencies. 

Note that the form of Pauli rate equation remains valid 
itself in the strong coupling limit. In the present paper 
we employ a perturbative approach to the rate equation 
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APPENDIX: RATE COEFFICIENTS 

Introducing the coupling reservoir spectrum density 

J a (u) = irJ2 q \ T g ~ e aq) and applying Fermi 

statistics to the reservoir modes, the correlation func- 
tions (0 and ((6l) can be written as 

C (±) (t) = r^j QM/ (±) He T^. (a.i) 

J-oo 7T 

Here, fL +) (tu) = 1 - fL~\uj) = [1 + e^"-^)]-! is the 
Fermi distribution function, with (3 = l/fc&T being the 
inverse temperature and /i a the chemical potential of to 
the a-electrode. Adopting a Lorentzian form of spectral 
density, 

W =-£/[(" -fi«) a +7a]> (A.2) 

the finite spectral width parameter ^ a is used to char- 
acterize the non-Markovian nature of system-reservoir 
coupling. With the complex roots of the Fermi function 
and of the Lorentzian spectral density, the integrals in 
Eq. (|A.1[) can be determined by the residues of the Ker- 
nel. The resulting infinite series are 

2? - 

ci +) (t) = j«fL +) (A a y At - -J2 Uv k y Vkt (a.3) 

p k=i 

and 

2 m 

CtHt) = j^i-A^e-^ _ _£j a („2) e -K«, 

p fc=i 

(A.4) 
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with the abbreviation A Q = £l a +i"f a and Vk = ^ a + i^k, 
where W}, = (2k — l)7r//3 are the Fermion Matsubara 
frequencies. In order to completely separate real and 
imaginary parts of the correlation functions (0 and ^ , 
we first separate its individual components. For the two 
complex Fermi functions we calculate 



1 + e ±/3(n a - Mc ,) e -i/3 7Q 



(A.5) 



where 

X± = l + 2cos(/3 7Q )e ±/3( ^-^ ) +e ±2^<»-n«). (A.6) 
The complex spectral densities can be separated into 



iw k ) 

[zu 2 k + (fiat - il a ) 2 + jI] =f 2im k {Ha - H a ) 



where 



— {^l + {^-^ a ) 2 +ll?+^l^c 



(A.7) 



(A. 



Based on the separation in real and imaginary contribu- 
tions, we can write the correlation functions as 

m 

fe=l 



The coefficients are all real: 



a n = 



X; 



:[l + e 



cos(/?7 Q )], 



X a 



^sin(/3 7a )], 



4(/i a — Q a )Tu k 

Wa ' 



(A.lOa) 
(A.lOb) 

(A.lOc) 

(A.lOd) 



Rigorously, the sum over the Matsubara values would be 
infinite; i.e., fc max — 77i — ► oo in Ecjs. (|A.3p and (|A.4p . 
but it can be truncated for practical purposes at a finite 
value that depends on the temperature of the system T 
and the spectral width. 
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